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In this work we study the spectral zeta function associated with the Laplace operator acting on 
scalar functions defined on a warped product of manifolds of the type I Xf N where / is an inter- 
val of the real line and N is a compact, cZ-dimensional Riemannian manifold either with or without 
boundary. Starting from an integral representation of the spectral zeta function, we find its analytic 
continuation by exploiting the WKB asymptotic expansion of the eigenfunctions of the Laplace op- 
erator on M for which a detailed analysis is presented. We apply the obtained results to the explicit 
computation of the zeta regularized functional determinant and the coefficients of the heat kernel 
asymptotic expansion. 



I. INTRODUCTION 



The spectral zeta function is one of the most widely used tools for the analysis of the spectrum of 
a class of partial differential operators generally defined on Riemannian manifolds. It is often the case, 
in both mathematical and physical problems, that particular information needs to be extracted from the 
spectrum of an elliptic self-adjoint differential operator with positive leading symbol defined on compact 
Riemannian manifolds with or without boundary. For the Laplace operator, in this situation, the spectrum, 
l n , is discrete, bounded from below and forms an increasing sequence of numbers tending to infinity with the 
behavior l n ~ w 2 / dim ( M ) where M denotes the manifold under consideration [22]. The spectral zeta function 
is constructed from the eigenvalues l n as 

oo 

fW = &'- (Li) 

11=1 

where s is a complex variable and each eigenvalue is counted with its (finite) multiplicity. Due to the 
asymptotic behavior of the eigenvalues mentioned above the representation of the spectral zeta function 
(11.11) is valid in the halfplane %(s) > dim(M)/2. It is possible, however, to analytically continue £(s) to a 
meromorphic function in the entire complex plane possessing only simple poles and which is holomorphic 



at the point s = [28]. 
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In physics the spectral zeta function is of pivotal importance because it provides an elegant way of 
regularizing the divergent quantities that often plague calculations performed in the ambit of quantum field 

nr 



theory in flat or curved spacetimes BIOL 



23h . Zeta function regularization techniques are predominantly used 



in order to compute, in a variety of situations, the one-loop effective action and the Casimir energy (see for 



instance (j, 



4DQI2I 



241P . For these types of applications one needs to evaluate the derivative of the 



spectral zeta function at s = or the value at s = - 1 /2 to obtain the functional determinant respectively the 
Casimir energy. Since these points do not belong to the region of convergence of (11.11) methods that provide 
the analytic continuation of (11.11) to values of %(s) < dim(M)/2 need to be developed. One of these methods 
relies on a complex integral representation of the spectral zeta function based on the Cauchy residue theorem 



II24I1 and it has been proven to be very useful in problems involving a wide range of geometries. It is this 
technique that we will employ here in order to find the analytic continuation of the spectral zeta function 
for Laplace operators on warped products of manifolds. 

Warped products of manifolds of the type Ix fN with / being an interval of the real line and N a compact, 
ci-dimensional Riemannian manifold play an important role especially in field theoretical models inspired 
by string theory. In the Randall-Sundrum two-brane model the spacetime is assumed to be five dimensional 
with only one 5 1 /Z2 orbifolded extra-dimension. The solution to the five dimensional Einstein equations 
which preserves Poincaire invariance is a warped product manifold with exponential warping function and 



four dimensional Minkowski branes 12911 . Generalized Randall-Sundrum models with higher-dimensional 



2 ill and references therein). In 



curved branes in an AdS bulk have also been considered (see e.g. HI, 
this paper we will mainly be concerned with warped product manifolds possessing an unspecified, strictly 
positive, warping function and a cf-dimensional smooth compact manifold N. 

The Laplace operator acting on scalar functions on a warped product IXf N is separable and therefore its 
spectral zeta function is suitably analyzed by using the contour integral method. Its analytic continuation is 
found by exploiting the asymptotic expansion of the eigenfunctions for large values of a specific parameter. 
For many geometric configurations considered in the literature the eigenfunctions are Bessel functions and 



their well known asymptotic properties have been used in order to perform such analytic continuation j2411 



Recently, it has been shown how to find the analytic continuation of the spectral zeta function in the ambit 
of the spherical suspension where the relevant eigenfunctions are associated Legendre functions [17]. The 
case of a warped product / x f N is more general because for an arbitrary strictly positive warping function 
the eigenfunctions are not explicitly known. In this paper we will show that even when the warping function 
and the manifold N have not been specified the asymptotic expansion of the eigenfunctions can be found 
and the process of analytic continuation can be carried to completion. This will allow us to obtain explicit 
formulas for the zeta regularized functional determinant and for the coefficients of the heat kernel asymptotic 
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expansion of the Laplacian on the warped product manifold. 

The outline of the paper is as follows. In the next section we describe the geometry of the warped 
product I Xf N and present the eigenvalue problem for the Laplacian. In sections [III] and [IV] we compute 
the asymptotic expansion of the relevant eigenfunctions and use it in order to explicitly perform the analytic 
continuation. The analysis needs modifications if the Laplacian on N has zero modes. These are described 
in section [V] Finally, we use the results of the analytic continuation in order to obtain the regularized 
functional determinant in section [VI] and the coefficients of the heat kernel asymptotic expansion in section 

eh 



II. THE SPECTRAL ZETA FUNCTION 



In this work we will consider a bounded D-dimensional warped manifold constructed as follows: Let 
/ = [a, b] c R. and let N be a (i-dimensional compact Riemannian manifold with or without boundary dN. 
For / e C°°(I), and f{r) > with r € /, we consider the warped product M = I Xf N by referring to / 
as the warping function. It is clear, from this construction, that M is a compact manifold with boundary 
dM - N a U Nb, where N x denotes the cross section of M at the point x € I. For the moment we will assume 
that / is a smooth function on the interval /, but depending on the type of information one wishes to extract 
from the spectral zeta function, this assumption can be relaxed. The local geometry of M is described by 
the line element 



As 1 = dr 2 + / 2 (r)dl 2 , 



(2.1) 



with dZ 2 being the line element on N. The operator we are interested in is the Laplacian acting on scalar 
functions tp e £}{M), namely 



AmV> = p — 9^ (f v Vdetg d v ) <fi . 
V det £ 



(2.2) 



By denoting with X a set of local coordinates on N, the metric tensor for the warped product M can be 
written in the form 



g(r,X) = 



1 

•2c 



(2.3) 



(0 f z (r)hij(X)) 

where hjj(X) is the metric tensor on the manifold Af and = {1, • • • , d}. By exploiting the metric (12.31) . 
the Laplace operator (12.21) can be written as 



A W<r ,x) = |-^ + ^A + _L_ AK| , (r , x) , 



(2.4) 
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where A# denotes the Laplace operator on the manifold N. 

In the following we will be interested in the pure Laplacian and we will therefore consider massless 

scalar fields minimally coupled to the curvature. In these circumstances the relevant eigenvalue problem 
reads 

- A M <f(r,X) = A 2 <p(r,X) . (2.5) 

An ansatz for a solution to the problem (12.51 ) can be found by separation of variables. Let "K(X) € £?{N) 
be the <i(v)-times degenerate harmonics on N satisfying the eigenvalue equation 

- Ai/H(X) = v 2< H{X) . (2.6) 

An eigenfunction of (12.41 ) can be written as a product 

<p(r,X) = tfrYHQO , (2.7) 

where (f>(r) satisfies the second order differential equation 

d 2 fir) d , v 2 \ 

^ +d -TTT^ +A \<Kr) = 0, (2.8) 



x dr 2 f(r) dr f\r), 
and the prime denotes, here and in the rest of the paper, differentiation with respect to the variable r. 
Since the warping function f(r) > for r e / and, in particular, /' € C (I), the coefficients of the above 
differential equation are continuous in / which ensures that a continuous solution cp(r) to ( 12.81 ) in the given 
interval exists. A unique set of eigenvalues and eigenfunctions of (I2.8I ). and hence of (12.5I ). is found once 
suitable boundary conditions are imposed. For defmiteness, we will consider both Dirichlet and Neuman 
boundary conditions which lead to the requirements 

cp{a) = cp{b) = , respectively (p\a) = cp\b) = . (2.9) 

The spectral zeta function for the Laplace operator on M is defined as 

A 

where we will assume that no negative eigenvalues occur so that we can, as is standard, use the non-positive 
real axis as the branch cut of the logarithm. The series that defines the zeta function in ( 12.101 ) is convergent 
for %(s) > D/2 and can be analytically continued, in a unique way, to a meromorphic function in the 
entire complex plane which coincides with (12.101 ) in its domain of convergence. It is the construction of the 
analytic continuation of ( 12.101 ) to values of %(s) < D/2 on which we will be focusing in the first part of this 
work. 
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Since N is unspecified, we will express the spectral zeta function on the whole manifold M in terms of 
the zeta function £n(s) associated with the operator -A#. Its definition is (9J] 

£ N (s) = Y,d(v)v~ 2s , (2.11) 

y 

which is convergent for %(s) > d/2. 

The analytic continuation of the spectral zeta function in (12.101 ) can be carried out by expressing its 
defining series in terms of a contour integral by exploiting Cauchy's residue theorem. Obviously, the eigen- 
values A are not known for a general warping function /. However, they can be determined in an implicit 
way. For this purpose, when Dirichlet boundary conditions are imposed, it is convenient to consider the 



following initial value problem jlo. 



25, 



26] 



dr z j(r) ar J \r)j H 

where p e C and u p (r, v) is an analytic functions of p. The eigenvalues A of the original boundary value 
problem 02.81 ) and 02.91 ) are recovered as solutions to the secular equation 

u p (b,v) = 0. (2.13) 

The condition 02.13b allows us to represent the spectral zeta function on M in terms of a contour integral 



which is well defined in the region %(s) > D/2, in more detail we have j5l47l. 1 13112411 

«j) = 2 rf Wv( f ). ( 2 - 14 ) 

V 

where 

Us) = 2^J c d P (P 2 + m2 Y Q-p ln Up{b ' V) ■ {2A5) 

In the above equation, C represents a contour in the complex plane that encircles all the roots (assumed to 
be positive) of u p (b, v) and we have also introduced, for technical convenience, a mass parameter m which 
is supposed to be small and will be sent to zero in the final results. 

By deforming the contour C to the imaginary axis and by using the property ui p (r, v) = u-i p (r, v) we 
obtain, once the change of variable p — » zv has been performed, the representation 

&(*) = — r°dz(vY-m 2 r ^\nu iyz {b,v) , (2.16) 

V 

which is now valid in the region 1/2 < %(s) < 1. It is important to point out that we work under the 
assumption that no zero modes of the Laplace operator -Ajv are present. If v - is instead an eigenvalue 
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of -An then the process of analytic continuation needs to be slightly modified as will be outlined in section 

El 

The case of Neuman boundary conditions can be analyzed by following a procedure similar to the one 
described above with the addition of a few modifications. For Neuman boundary conditions it is necessary 
to replace the problem in (12.121 ) with 

^2 f'(r) d 2 v 2 \ 

"7~7 + d ~T7T j — *~ P 2 — 7T7T ) u p( r > v ) = , u p (a, v) = 1, u' {a, v) = . (2.17) 
dr 2 f(r) or f in) 

The original eigenvalues A of (12.81 ) are then implicitly obtained from the equation 

u' p (b,v) = 0. (2.18) 

The process of analytic continuation follows the same lines as described above leading to the expression for 
the spectral zeta function in the Neuman case 

£ N (s) = J]d(v)tf(s), (2.19) 

V 

where 

^v%)-— J m d Z (vV-m 2 ) -\nu' ivz (b,v) , (2.20) 

V 

which is well defined in the region 1/2 < %(s) < 1. 

In order to obtain the analytic continuation of the integral representation (12.161 ) to values of ^R(s) < D/2 
we add and subtract a suitable number of terms of the uniform asymptotic expansion of the functions 
Ui vz {b,v) for v — > oo and z - p/v fixed. The same procedure of analytic continuation can be applied 
to (12.201 ) by utilizing, instead, the uniform asymptotic expansion of the derivative u' iv (b, v) for v — > oo and 
z = plv fixed. In the next section we will provide explicit expressions for the uniform asymptotic expansions 
needed for our analysis. 



III. THE WKB APPROXIMATION OF u ivz AND u\ vz 

The differential equation (12.121 ) can be transformed into an equivalent one which does not contain the 
first derivative term. By defining 



U(r) - J^TT ' and y ( y > z ' r ) = ~ y2 ( zl + 727T I » ^ 
f(f) \ f 2 (r) 



the above mentioned transformation is achieved with the use of the ansatz 

u ivz (r, v) - exp j-i JT U(t) df J ¥ v (z, r) . (3.2) 
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By substituting (13.21 ) into (12.121) . with p 2 = -v 2 z 2 , we obtain the following differential equation satisfied by 
¥ v (z,r) 

u 2 \ 

— I + q(v,z,r)W v (z,r) = 0, (3.3) 
where we have introduced the function 



q{v,z,r) = V{v,z,r)-^U\r)-^U 2 (r) 



-v 2 z 2 + 



1 



df"(r) d{d-2)f'\r) 



f 2 (r)) 2 fir) 



4 f\r) 



(3.4) 



The large-v asymptotic expansion of ( 



3.3b. valid uniformly with respect to the parameter z, can be ob- 



tained by exploiting the WKB technique (2j,|27j]. It is well known that the uniform asymptotic expansion of 
the above equation for v — > oo contains exponentially growing and exponentially decaying terms. Although 
for the purpose of the analytic continuation of the spectral zeta function we will only need the exponentially 
growing part, at the beginning we have to consider both types of contributions in order to be able to impose 
the initial condition in (12.121 ) or in ( I2.17I ). The introduction of the auxiliary function 



S(y,z,r) = — ln¥ v (z,r) 
or 



(3.5) 



leads to the non-linear differential equation 



S'(v, z, r) = -q(v, z, r) - S 2 (v, z, r) , 



(3.6) 



which is obtained by differentiating S(y,z, r) and by noticing that *¥ v (z, r) satisfies (I3.3I ). For v — > oo we 
seek an asymptotic expansion for S(v, z, r) of the form 

" Si(z,r) 



S(v,z,r) ~ vS-i(z,r) + S (z,r) + ^ 



i=i 



v' 



(3.7) 



which is valid uniformly with respect to the parameter z. By substituting the above ansatz into (13.61 ) and by 
equating like powers of v we obtain explicit expressions for the terms S i(z, r) with i = {-1, 0, 1, • • • }. For 
the leading and the first subleading term of the asymptotic expansion (13.71) one finds 



St,(z,r) = ±1 2 + , S± (z,r) = -lllnSt^r) 



(3.8) 



For the terms proportional to v ' one has, for i = 1, 



Sj(z,r) = - — 



1 



df"(r) d(d-2)f'\r) _ c2 



2 f(r) 4 f\r) 



+ S l (z,r) + S' Q (z,r) 



(3.9) 
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and the following recurrence relation valid for the higher asymptotic orders with i > 1 

1 



Sf + Az,r) = - 



25 f, for) 



5 ^fo r) + $ « fo r)S f_ n (z, r) 



n=0 



(3.10) 



The ± stands for the choice of sign in the leading order term of the asymptotic expansion, namely S-\(z, r), 
and it will give rise to two independent solutions S + (v, z, r) and S~(y, z, r) to the differential equation ( 13.61 ). 
These solutions, thanks to the relation (13.51) . will determine the functions *F+ fo r) and r) which, in 
turn, will give the exponentially growing and decaying behavior of the functions w,v z (r, y) (cf. (13.21) ). 

The correct large-v asymptotic expansion of the solution ( 13.21 ) is given in terms of the linear combination 



Ui vz {r, v) = exp 



_i r 

2J a 



1/(0 df 



J a 



Aexpl S + {v,z,t)dt \ + Bexpl S~(v,z,f)dt 



J a 



(3.11) 



where the constants A and B are uniquely determined once the initial conditions are imposed. 



A. Dirichlet Boundary Conditions 



The initial conditions associated with the Dirichlet case, namely Uj vz (a, v) = and u' t (a, v) = 1, imply 

1 



A +B = , A = 

o + (v, z, a) - S (v, z, a) 

By substituting the explicit expression (13.121 ) for A and B into (13.1 II ) we obtain 



(3.12) 



Uiv Z (r, v) 



— -i— -expj-i f U(t)dt\exp{[ S + {v,z,t)dt\(l+B{v,z,r)), (3.13) 



where the function £(v, z, r) contains the exponentially decaying terms of the asymptotic expansion. For the 
purpose of the analytic continuation of the spectral zeta functions (12.141 ) we actually need the asymptotic 
expansion of In u ivz (r, v) which from 03.131) follows as 

In u ivz (r, v) = - In [S + (v, z, a) - S~(v, z, a)] - )- f 17(0 dt + f S + (v, z, t)dt + £(v, z, r) , (3.14) 

^ Ja Ja 

where £>(v,z, r) denotes exponentially small terms as v — » oo. 

The formula obtained above, together with the results ( 13.71 ) through ( 13.101 ), will allow us to find explicitly 
the uniform asymptotic expansion of In w,- vz (r, v). Let us start with the analysis of the first term in (13.141 ). 
From (I3.8I) - (I3.10I) it is not very difficult to obtain the expansion 



In [S + (v, z, a) - S~(v, z, a)] = In (2v) + i In ( z 2 



1 \ °° 



O/foa) 



j+i 



(3.15) 



where the a) are determined through the cumulant expansion 



In 



, l/ 2 I " ~ 
1 + r r + 



E 



0Jk(z, a) 



k+\ 



Offoa) 



v i+\ 



(3.16) 
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and we have introduced the functions 



a>i(z,a) = SUz,a) - 5, (z,a) . 



(3.17) 



The result (13.151 ) and the expansion (13.71) corresponding to S + (v, z, a) provide an explicit expression for the 
uniform asymptotic expansion of In Ui vz (b, v), namely 



In Ui VZ (b, v) - - In (2v) - - In I z L + 



1 



1 



. H — 7 In 

/2(a)/ 4 



1 + z l f l {a) 
1 + Z 2 f 2 (b) 



d-\. m 
+ — ln W) 



+ v 



Ja tl V 



(3.18) 



where we have discarded exponentially small terms and we have exploited 



f + If 11+ z 2 f 2 (a) 

In addition, the functions jVt;(z, a, &) introduced in (13.181 ) are defined to be 

M l {z,a,b)= f S\{z,t)dt, 

Ja 



d-l f(a) 
+ In 

2 f(b) 



(3.19) 



(3.20) 



while for i > 2 



Mi(z, a, b) 



Jo 



t)df-Di_i(z,a) . 



(3.21) 



Before proceeding with the study of the analytic continuation of the spectral zeta function we need to find 
an expression for the functions Mi(z, a, b) which shows in an explicit way their dependence on the variable 
z. This is necessary in order to be able to perform the z-integration present in (12. 16b - From the recurrence 
relation (13.101) and from the analysis of the first few St(z, r) one can deduce the following general form 

i+i 



F k dr) 



k=0 { V 2 + _L- 

P(r) 



2k+i 
2 



(3.22) 



The dependence on the variable z is now manifest, however a relation for the newly introduced functions 
Fk, i( r ) is needed. From the explicit expression for Si (z, r) and S J" (z, r) in (13.81 ) and (13.91 ) one readily obtains 

f'(r) 



*o,oto = 0, F l0 (r) = 



2f\r) 



(3.23) 



df'ir) d(d-2)f\r) 
Fo,i(r) = -r—r— + 



4 m 



3f'\r) If "(f) 



f 2 (r) 



4 f\r) 4 f\r) ' 



F 2 ,i(r) = -- 



5/' 2 (r) 



8 f6( r ) 



(3.24) 



(3.25) 
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The functions Fkj(r) of higher order are found by substituting (13.221 ) into the recurrence relation (13.101 ) and 
by equating like powers of (z 2 + f~ 2 (r)), in more detail for i > 1 



Fk,i+i(r) = -- 



(2k + i - 2) ^^ (r) + F' k ( .(r)#(i + 1 - k) + n, i(r) 
P(r) 



where H(x) denotes the Heaviside step function and ,(r) is defined through the equality 



i m+li-m+1 , ^ 2*+2j+/ ;+2 

m=0 it=0 ;=0 \ J \ )l n=0 



Z + 



1 



f\r) 



(3.26) 



(3.27) 



It is important to stress that the relations (13.261 ) and (13.271 ) can be easily implemented in an algebraic com- 
puter program in a way that allows one to compute Fkj(r) to any desired order. For convenience the 
functions Fkj(r) are given in appendix lAl up to the order i = 3. 

Let us turn our attention to the functions D;(z, a). From the recurrence relation (13.101) one can show that 
the following holds for i > 1 



S;(z,r) = (-\ySf(z,r), 



which implies, in particular, that 



a>2i(z, a) = , and oj 2 i+\ (z, a) = 25 J +1 (z, a) . 



(3.28) 



(3.29) 



The above expression shows that only even inverse powers of v will appear in the expansion (13- 16b - Inspec- 
tion of 03.291) and the cumulant expansion ( 13.161 ) suggests the following general form for i > 1 

Q. k j(a) 



f(a) 



k+i ' 



(3.30) 



where the functions ,(a) can be found from the relation 



In 



oo j ' 2j i . - n - j\ co ^ 2i I ^ \-k-i 

'^Sil^WP^ "^S**T + ?M (33!) 

j=l \n=0 x J ' )\ i=l k=0 x J " 



by equating like powers of v and (z 2 + f~ 2 (a)). We have finally arrived to the point in which we can write 
an expression for the functions Mi(z, a, b) that explicitly shows the z-dependence. Since Dij(z, a) = with 
j € N + we need to distinguish between two cases: when i = 2m + 1 with m € No, then 



M 2m +i(z,a,b) = 
while for i = 2m with m e N + we obtain 

2m+l n b 



2m+2 n \, 

^ / 

h-(\ J a 



Mim(z, a, b) 



dtF kt 2m(t)\zf + 



dtFk,2m + l(t)\Z 1 + 



1 



1 



fHt) 



fHt) 



-k-m 2m 



- ^jQk,m(a) 



k=0 



z 2 + 



1 



f 2 (a) 



-k-m 



(3.32) 



(3.33) 



The last two formulas, together with the expansion ( 13.181 ), represent all the information we need in order to 
perform the analytic continuation of the spectral zeta function. 
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B. Neuman Boundary Conditions 

For Neuman boundary conditions the relations u p (a, v) = 1, u'Ja, v) = lead to the following solution 
for A and B in (f37TTb 



A = 
B = 



1 



\u(a)-S~{v,z,a)\ , 



<S + (v, z, a) - S (v, z, a) \2 

— -i— -IS + (v, Z ,a)-^(«)| , 

<M(y,z,flD - <i> (y,z,a)\ 2 



(3.34) 



By substituting the above formulas into the ansatz ( 13.111 ) and by differentiating the resulting expression we 
obtain the uniform asymptotic expansion 

[S-(v,z,a) - \U{a)\ [S + (v,z,r) - \U{r)\ 



<S + (y,z, a) - S (v,z,a) 



x exp U(i) dfj exp jj" S + (v, z, f)df j (l + ft(v, z, r)) , (3.35) 

with ^(v, z, r) denoting exponentially small terms. By comparing (13.351) with (13 - 1 3b . it is straightforward to 
derive the expression 



In u'j (b, v) = In Ui vz (b, v) + In 



-S-(v,z,a) + -U(a) 



+ ln 



S\v,z,b)-\u{b) 



(3.36) 



where the exponentially small terms have been neglected and we have set r - b. From the expansion (13.7 
and the results (13- 8b - (l3- 10b we have, for the last two terms in (13.36I ). 



In 



-S-(v,z,a) + -U(a) 



1 , 

In v + - In z + 



/ 2 (a)) + § 



Z;(z„a) 



and 



In 



S + (v,z,b)-lu(b) 



= In v + - In I z z + 



f\b) 



(=1 



(3.37) 



(3.38) 



The functions "Z,. (z, a) and Z^fo b) are defined through the cumulant expansions 



In 



1+ z 2 + 



f 2 (a) 



Y ^ a ) 



;=1 



and 



In 



1+ z 2 + 



-fx 



f 2 (b) 



<^,b) 



k=\ 



i=l 



(3.39) 



(3.40) 



with 



o- Q (z, a) = -U(a)- S (z, a) , cr k (z, a) = -*S A (z, a) , 
cr+(z, ft) = <S+(z, fe) - i U(b) , <x+(z, b) = S + k (z, b) , 



(3.41) 
(3.42) 
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where k > 1 . 

The functions cr~(z, a) and crt(z, ft) can be written explicitly in terms of the variable z by comparing their 
defining formulas in (13.411 ) and (13.421) with (13.221 ) and (13.28I ). In more detail we have 



cr 



"fc «) - E 2177 . <fe « = Ya 



k=0 | 2 + 1 



- ( z2 + m) 



2k+i 
2 



where, with i > 1 and < k < i + 1, 



and 



(3.43) 



(3.44) 



K o(b) = -■= U(b) , <U\ (b) = Fi, (6) , <W+ ,.(ft) = F*, ; (ft) • 



(3.45) 



From the expression (13.431) and the expansion ( 13.391 ) one can show that the functions Z. t (z, a) and Z*(z, b) 
have the form 



z ; (z,a) = 2j mr - z i (z,b) = 2_ J 



k=Q 



k=Q 



(«- 2 + ?k) V 



where <3, (a) and <3^ .(b) can be found through the relations 



In 



2k+i 
2 



(3.46) 



(3.47) 



In 



yJ 

j=l \ n =Q 



1 \ — ^ 



Pib) 



f 2 (b) 



2k+i 
2 



(3.48) 



by equating like powers of v and (z 2 + /~ 2 ). 

By substituting ( 13.391 ) and (13.401 ) into the relation (13.361 ) and by making use of the uniform asymptotic 
expansion (13.181 ) we obtain 

l+z 2 / 2 (a)" 



\nu' ivz (b,v) = 21nv-ln(2v) + iln|z 2 + ^J + iln 



+ v 



S^iz, t)dt + V - [M(Z, ^ + ") + ^f"fe ^] . 



rf-l, /<«> 

+ — ln 7w 



(3.49) 



which will be used for the analytic continuation of the spectral zeta function when Neuman boundary 
conditions are imposed. 
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IV. ANALYTIC CONTINUATION 



The starting point of the analytic continuation for the spectral zeta function associated with Dirichlet 
boundary conditions is the integral representation (12.161) . By adding and subtracting L leading terms of the 
uniform asymptotic expansion of the functions found in (13.181) we obtain 

I 



(is) = Z(s) + £ A f (5) , 



t=-l 



where 



Z(s) 



sin jts \ i r°° / 9 9 9 \-.v d { , 1 / 9 1 

— 2 d(v) J a dz (v 2 z 2 -m 2 ) -| In uUb, v) + In (2v) + - In | Z 2 + — 



Iln 



1 + z 2 / 2 (&) 
and the functions A,(j) are 



d-1 /(a) 



In : — — v 

/(*) 



M,(z, a, fc) 



A. 1 (s) = ——2_ J d(v)J m dz(v 2 z 2 -m 2 ) - vj S^feOd* 



(4.1) 



(4.2) 



(4.3) 



z 2 + ^) + Iln 

/ 2 («)/ 4 



(4.4) 



Mi(z, a, b) 



(4.5) 



The function Z{s) in (14.21 ) is analytic in the region %(s) > (d - 1 - L)/2. This restriction is found by 
considering the behavior on the integrand in v, which is v - 2s - L - 1 ^ and by recalling that the spectral zeta 
function on N converges for %(s) > d/2 according to Weyl's estimate. The integrals appearing in ( 14.31 )- 
(14.51) are well defined, for any fixed v, at least in the strip 1/2 < %(s) < 1. At this point of the analysis 
the assumption of smoothness of the function f(r) can be relaxed. Let / e C L+2 (I) and /(f) > for 
tel. Under this assumption the integral of S^^zJ) over the interval / is uniformly convergent for z e 
[m/v, oo), and in addition the derivative d z S^(z, t) is continuous. This implies that we can interchange in 
(14.31) the differentiation with respect to z and the integration over the interval /. Moreover, since d z S^{z, t) 
is continuous for (z, t) € [m/v, oo) x / we can interchange the integration over the variable z with the one 
over the variable t for the values 1/2 < %(s) < 1. Furthermore, since the functions S?(z, t) contain at most 
the derivative of order L + 1 of the function / the same remarks outlined above apply also, for 1 < i < L, to 
the integrals in ( 14.51 ). 
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By differentiating 5^(z, t) with respect to z and by performing the change of variables z — > m/(uv) we 
obtain 

A y is) = S ^l m -2s + l r du u 2s-2 (l _ u 2ys f" & L + . (4.6) 

«" Jo J« \ m 2 f 2 (t)l 

The integration over the variable u, for 1/2 < %(s) < 1, leads to the result 



i^u k\ m I 2jj a 



By using the binomial expansion for mf(t)/v <sc 1 and by subsequently summing over the angular eigenval- 
ues v we obtain the final expression 

i f(-i) tr ('+*-i)^ ar J. 1 i\ f* 

where we have used the definition ( 12.1 11) . By following the procedure just described one can find expressions 
similar to (14.81 ) also for Aq(s) and A;(j). In particular, for Aq(j) one has 

A ^ = -j X ^^rrr™ 2 ^ + *> [/ 2,+2 *^ + / 2,+2 *^l • ^ 

In addition, for i = 2n + 1 with neNo 



a 2 „ + i(^) = 2j — ^ r + + " + 2 j 

2 p?r(s + j + n + k+ i) r fc 

E-^-T n^M ^,2„ + i(0/ M+2i+2 " +1 (0d ? , (4.10) 



while for / - 2n with n£N + 

1 °° (-l) k 
Amis) ■ L -ry-rn 2k ^ N {s + k + n) 

2/1+1 



% Yijll) \l F jMt)f S+2k+2j+2n m - £1 7 ; n (a)f 2s+2k+2J+2n (a) j , (4.11) 
where fi2n+i,n(«) = 0. 

The expression (14.11) with the results (14.81) through (14.111 ) represent the analytic continuation of the spec- 
tral zeta function f (s) on the manifold M and it is written in terms of the zeta function, £/v(s), on the manifold 
N. The function Z(s) is analytic in the region %(s) > (d- 1 -L)/2 while A ; (j) are meromorphic functions in 
the entire complex plane. The choice of L, namely the number of leading terms of the asymptotic expansion 
of In Ui vz (b, y), depends on the specific quantity that one wishes to compute by using g(s) [24]. In particular, 



in order to evaluate the functional determinant and the anomalous scaling factor of the operator -Am one 
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needs the value £'(0) respectively £(0). This means that it is sufficient to add and subtract L = d leading 
terms of the asymptotic expansion. For the calculation of the «-th coefficient of the heat kernel asymptotic 
expansion for Am it is sufficient, instead, to set L = n - 1. 

For Neuman boundary conditions the analytic continuation proceeds along the same lines described 
above. From the integral representation (12.201 ) and the uniform asymptotic expansion (13.491 ) one obtains 

I 



£ N (s)=Z N (s)+Y J A?(s), (4.12) 
i=-i 

where Z N (s) is analytic in the region %(s) > (d - 1 - L)/2 and has the representation 

Z"( S ) = ^^J(v) j^°°dz(yV-m 2 )" i ||ln M ; v ^,y)-21nv + ln(2v) 



Uh + ' -Iln 



2 \ P(b)! 4 

L 



l+z 2 f 2 (b) 



d — 1 f(a) , 
—— ln^-v S\(.z,f)dt 



fib) 



J a 



1 1 

" Z ~i [ M i( z ' a >® + + Z!(z,b)] | . (4.13) 

(=1 ' 

In addition, by comparing the uniform asymptotic expansions (13.181 ) and ( 13.491 ) one can show that the 
functions A. (s) in (16.11) can be written in terms of Aj(s) as follows 

A^Cs) - A_!(j) , A^(j) = -A Q (s) (4.14) 

and, for i > 1, 

A ^ s) = A ^-W)t ^ m M s+k+ i) 

2 F ^ + ^ + ^ + ^ [Q • ; («)/ 2i+2 ' +2y+, '(«) + flj ^1^(4)] • (4- 15) 

7=0 T \j + 2) 



X 



The results (14.131 ) through (14.151 ) provide the analytic continuation for the spectral zeta function for Neuman 
boundary conditions. 



V. PRESENCE OF ZERO MODES ON THE BASE 



The analytic continuation of the spectral zeta function illustrated in the previous section was performed 
under the assumption that no zero modes on the base were present. In this section we will consider the 
case in which v = is an eigenvalue of -Aa? with degeneracy d(0). In this situation the process of analytic 
continuation needs to be amended in order to take into account these modes. The integral representation 
of the spectral zeta function for Dirichlet boundary conditions is obtained by separating the contribution 
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coming from v = 0. Explicitly it reads 

m =^~Tj d{v) x dp ( p2 + mi y Tp in up{h ' v) + 1*1 f c dp ( p2 + m2 y s t p in kp(b) ' (5- 1} 

where the sum is over v > and the function h p (r) is the unique solution of the initial value problem 
d 2 J'(r) d 2 



d; .2 + d J^d? + P J = ' ^ = ' ^ (fl) = 1 ■ (5 - 2) 

The analytic continuation of the first integral in d5.ll ) has been developed in the previous sections and will 
not be repeated here. We will, therefore, only focus on the analytic continuation of the second integral in 
(TO) . 

By deforming the contour of integration to the imaginary axis we obtain the following representation 

sin7rs T 1 / 1 9 \-s d smns C°° / 9 9 \-.v d 

&(s) = d(0) dz(z 2 -m 2 ) —lnh iz (b) + d(0) d Z (z 2 -m 2 ) — In h k (b) . (5.3) 

X J m v ' oz n Jy v 1 dz 

The first integral in (15.31 ) is valid in the region ^R(s) < 1 and the second integral in the region *R(s) > 1 /2. 

In order to be able to extend the validity of the above representation to the left of ^R(s) = 1 /2, we will need 

to add and subtract, from the second integral, the asymptotic expansion of the function h K {r) for z — > °°. 

The desired asymptotic expansion is obtained along the same lines of section |III] The function /z !z (r) is a 

solution of the problem (15.21 ) with p — > iz and by using the ansatz 

fc fe (r) = expj~i J [/(f)d/Ufcr), (5.4) 

we obtain the auxiliary equation 

d 2 



j2 +e(v,z,r)|O(z,r) = 0, (5.5) 



where we have introduced the function 



n< ^ 2 df"{r) d(d-2) f' 2 (r) 

Q(z,r) = -z - -— -pr— . (5.6) 

2 f(r) 4 f\r) 



Using the notation 



P(z,r) = ^-\n^(z,r), (5.7) 
or 

we have the relation 

P\z, r) - -Q(z, r) - <P 2 {z, r) , (5.8) 
and we seek an asymptotic expansion for z — > °o of the form 

°° P (r) 

P(z, r) ~ Zj P_!(r) + P (r) + J] ■ (5.9) 



;=1 
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By substituting (15.91 ) into (15.81 ) and by equating like powers of z we obtain 



P±(r) = ±l, P Q (r)=0, Pf(r) = ±- 



df"(r) _ </(</- 2) /' 2 (r) 



4 /(r) 



/ 2 (r) 



and the recurrence relation for i > 1 reads 



m=0 



(5.10) 



(5.11) 



By using ( 15.91 ) and the relation (15.71 ) into (15.41 ) we obtain, after imposing the initial conditions in (15.21) . 
the asymptotic expansion of the logarithm of the solution hj z {b) as 

In h iz (b) = -ln[P + (z,a)-'p-{z,a)] - I f tf(f)df+ f ^ + (z,0df, (5.12) 

^ Ja Ja 

where exponentially small terms have been omitted. From the asymptotic expansion (15.9I ). the results (15.101 ) 
and (15.111 ). and by noticing that Pj(r) = (-l)'Pt(r), it is not very difficult to show that we have 



In [P + {z, a) - P~(z, a)] = In 2z + ^ , 



(5.13) 



where the terms 7~j(a) are obtained from the cumulant expansion 



In 



1 + 



2 



Pu-\(a) 



2 



(5.14) 



We can finally write down an explicit expression for the asymptotic expansion of In hj z (b) valid when z — > oo. 
By exploiting (15.131 ) and ( 15.91 ) through (15.111 ) we obtain 



d f(a) \— 1 
In h iz (b) = -ln2z+-ln J j-±+z(b-a) + 2_ J 



2 f(b) 



where 



and, when i > 2, 



Ni(a,b) 



(0d/ 



(5.15) 



(5.16) 



(5.17) 



The analytic continuation of £q(s) in (15.31 ) is obtained in the same way as explained in section [IVJ and 
once the elementary integration in z of the leading asymptotic terms is performed we obtain, in the massless 
case, 



£o(s) = d(0) 

L 



sin ns P° 
n Jo 



dzz~ 2s — i \nh iz (b) + H(z- 1) 



d 



-2 



Ni(a,b) 



i=l 



Z< 



- d(0)- 



dz 

sin ns 



d f(a) 
In 2z - - In — — - z(b - a) 
2 fib) 



71 



2s 2s 



-T + V— 

1 4—* 2^ + i 

i=i 



(5.18) 



18 



where the integral represents an analytic function in the region %(s) > -{L + l)/2. 

When Neuman boundary conditions are imposed we obtain the following relation for the derivative of 
the eigenfunctions (cf. 13.361 ) 



InhUb) = In h iz (b) + In 



-p-(z,a) + ^U(a) 



+ ln 



9\z,b)- l -U{V) 



(5.19) 



The last two terms in the previous equation can be cast in the form 



In 



-p-(z,a) + -U(a) 



+ ln 



>P + ( z ,b)--U(b) 



CO 

21n Z + ^ - [Dj(a) + £>t(fc)] , (5.20) 

i=i z 



where the functions D- (a) and Df(b) can be found through the relations 



In 



I + Uifl) _ y 1 i w ^ y 
h 2z 4-* z i+l f-f z 



Pjia) 



D~(a) 



In 



Ujb) y * j \v> _ y 



2z 



i=l 



<DUb) 



(5.21) 



k=\ 



By utilizing (15.151) and the expansion (15.201) in ( 15.191 ) we obtain 

d f( ) °° 1 

\nh' iz {b) = 2\nz-\n2z+-\n J -^+z{b-a) + Y J - [^iia,b) + Djia) + £>+(£)] 



(5.22) 



The contribution of the zero modes on the base Af to the spectral zeta function for the Neuman case is 
obtained in the same way as for the Dirichlet case and it reads, for m = 0, 

.8 



d 3? ln/ 4W + #(z- 1) 



Sin 7TS C 

$(s) = d(0) 

n Jo 

L 1 11 

- zib - a) - ^ -. [Niia, b) + D7(a) + £>+(&)] I 

;=i z 



d f(a) 
In 2z - 2 In z In 

2 f(b) 



+ d(0)- 



smns 



(=1 



(5.23) 



where the integral is analytic in the region %is) > -(L + l)/2. The results (15.181 ) and (15.231 ) correspond to 
the additional contributions to the spectral zeta function that need to be taken into account when zero modes 
on the base are present and Dirichlet respectively Neuman boundary conditions are imposed. The important 
feature of expressions (15.181 ) and (15.231 ) lies in the fact that they render manifest the meromorphic structure 
of £ois) and $ (j). As a consequence, it is straightforward to extract, for instance, the residues which are 
the relevant quantities for the computation of the coefficients of the asymptotic expansion of the heat kernel. 



VI. ZETA REGULARIZED FUNCTIONAL DETERMINANT 



In this section we will utilize the results obtained regarding the analytic continuation of the spectral zeta 
function in order to compute the functional determinant of the operator - Am- We will be mainly concerned 
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with the massless case, however small mass corrections can be computed from the expressions (14.21 ) and 
(I4.8I) - (I4.11I) by following the techniques described, for instance, in [15]. In the framework of zeta function 
regularization the functional determinant of an elliptic, self-adjoint operator on smooth compact manifolds 
is defined as [ 10, 23, 30] 



det(-A M ) - exp{-<T(0)} , (6.1) 

where the derivative at the point s = is intended to be taken after performing a suitable analytic continu- 
ation of the spectral zeta function. It is clear, from (16.11) . that for the analysis of the regularized functional 
determinant we only need to compute £'(0). 

Let us start with the Dirichlet case and consider the function Z(s) in (14.21) . By choosing L - d, Z(s) is 
an analytic function for %(s) > -1/2 which implies that its derivative at s = is simply 

Z'(0) = - J^(v)jm« im (6,v)-v^ S! 1 (f,^)df + ln2v + iln|^ + -^ 

*Mi(f,a,bh 



d-l, f(a) 1, 

— ln M~4 ln 



1 + ^f\a) 



1 + 4f 2 (b) 



,2\ 



By expanding for small m and disregarding terms of order 0(m ) we obtain the following expression valid 
in the massless case 

Z'(0) = -^(v)jln Mo (£,v)-v£ r\t)At + \nlv-^-\nf(a) + ^-\nf(b) 

-X^'";"-"). (6.3) 

(=1 v > 

This represents the most explicit formula that one can obtain without specifying the warping function f(r). 
Moreover, even when a choice for f(r) is made, the above expression can only be treated numerically unless 
the warping function belongs to a special class for which the eigenfunctions of ( 12.121 ) are explicitly known. 

Next, we will focus our attention to the computation of the first derivative of the functions Ai(s) in (14.81 ) 
through (14. 111 ). For this purpose it is important to display the analytic structure of the spectral zeta function 
£)v(s + a) in the neighborhood of s - 0. It is well known that as s — > the zeta function ^{s + a) possesses 
simple poles for a — (d - k)/2 with k - {0, • • • ,d — 1} and for a = -(21 + l)/2 with I e No- This implies 
that when a coincides with a singular point we have the following Laurent expansion as s —> 



£ N (s + a) = -Res^V(a) + FP^(ar) + O(s) , + a) = ~-Res£ N (a) + O(s ) , (6.4) 

while when a is a regular point we have 

Ms + a) = t N (a) + s£ N (a) + 0(s 2 ) , (6.5) 
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where Res denotes the residue of the function and FP its finite part. By performing the derivative with 
respect to s and by taking into account the behaviors (16.41 ) and (16.51 ) of £n(s) we obtain at s = 0, in the 
massless case, 



A' (0) 



4(ln2-l)Res^V|-i|-FP^(-i 



J a 



(t)dt - 2Res £ 



N 



In f(t) 
fit) 



At , (6.6) 



and 



A o(°) = -\^ N (0) - ^jv(O) [In /(a) + In f(b)] . 



(6.7) 



The explicit results for the derivative at s = of the remaining Aj(s) are as follows: for i = In + 1 with 
n e No we have 

2h+2 T / , \ ( j+n , \ 



4n + i(0) = J] 2Res^ n+- In2-J — 
y=o L ^ ' v *=o ; 

Xfc / 1 \ 2;i+2 „fo 

^ 2n + i(t)f 2j+2n+l (t)dt - 2Res & ^ + - J ^ J Fj, 2 „ + i(t)f 2j+2n+l (t) In /(*)df 



(6.8) 



and when i = 2n with n eM + one finds 

2n+l 



)/ 2;+2 » 



A 2„(°) = - Z [ FP & + RCS & Vl] F J. 2n(t)f j+2n (t)dt - Qj, n (a 

j=Q { J « 

2n+i , „b \ 

- 2Res & («) ^ { J F i ^{t)f i+2n {i) In f(t)dt - Q,, „(a)/ 2 ' +2 » In /(a) j , 



(6.9) 



with H n denoting the «-th harmonic number. A formula for £'(0) which, in turn, gives the functional 
determinant for the Laplace operator Am through (16.11 ) can be easily obtained thanks to the relation 



r(0) = Z'(0) + A^(0) + A (0) + £ A^. +1 (0) + £ A^.(0) , 



(6.10) 



/=() 



and it reads 



f(0) = Z'(0) + Resf K |-- 



fit) 



dt 



~FP£v 



{t)dt 



11 12 2 { I 1\ 

- ~ 4(0) - 5^(0) [i n /(«) + ln f(b)] -ZZr ( 4' + ?r j; 2,+1 (a ' b) 

-Resfri + - 21n2-g — 



Qj,2i+\(a,b) - Ij,2i+i(a,b) 



[i] 2i+i ^ x 
J] 2 FP ^ (j) Q h 2i (a, b) + Res & (i) 2 , + i (a, *) - J/, 2/+i(fl, 6)] | , 

i=l 7=0 ^ > 



(6.11) 
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where for typographical convenience we have defined the functions 



Qjj{a,b)= FjjVfJ+XWt 

J a 



1 + (-!)' 



H ;; ,{a)f 2 i +i {a) , 



(6.12) 



and 



Ij, i(a, b ) = 2 \ a F j, i(t)f 2j+ \t) In f(t)dt - 2 1 + i 2 ly Q. , (a)fj +i (a) In f(a) . (6. 13) 



Another quantity of particular interest especially in quantum field theory is the anomalous scaling factor 
which, for a conformally invariant field propagating on a curved spacetime, coincides with the conformal 
anomaly. This term comes from the measure dependent part of the one-loop effective action and it is 



proportional to £(0) (see e.g. 1130). The analytic continuation obtained in section [TV] allows us to compute 
£(0) in an analytic way in terms of the zeta function on the manifold N. From the expression (14.21 ) by 
choosing L = d one can immediately set s = and find that Z(0) = 0, which means that there is no 
contribution coming from ( 14.21 ) to the value of the zeta function at s = 0. The only contributions to the 
anomalous scaling factor come from A,-(0) with -1 < i < d. In fact, from the explicit expressions (16.61 ) 
through (16.91 ), with m = 0, and by taking into account the expansions (16.41 ) and (16.51) we obtain 

/ 1\ C h 1 d I i\ M 

m = -Res £v J /" l (t)dt - -&(0) - g Res & (- j g Q h t {a, b) . (6. 14) 

Let us focus our attention to the case of Neuman boundary conditions. By setting L - d in ( 14.131 ) the 
function Z N (s) is analytic for ^K.(s) > -1/2 and its derivative at s = reads, in the massless case, 



(Z N )'(0) = - ^ d(y)i In u' (b, v) - v f r\t)At + ln2v - 2 In 
d 1 1 

1-1 ' 



J - 1 rf+ 1 

v-— — ln/( fl ) + — — ln/(fe) 



(6.15) 



The derivative of the asymptotic terms A J r(s) can be obtained in terms of the derivative of A,(j) using the 
relations (14.141 ) and (14.151) . More explicitly we find 



(A^)' (0) = A'_!(0) , (A^)' (0) - -A' Q (0) 



(6.16) 



and for / > 1 
(Af)'(0) 



= 4(0) - 2 {Res £v ( [r + y{j + 5)] + FPfr (ij) [QT.^f^ia) + G^)/ 2 ^)] 
;'=o 

- 2Res Cn (0 j • ,(«)/ 27+, '(«) In /(a) + <2+ i{b)f^\b) In /(&)] . (6. 17) 
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By differentiating (16.11 ) and by exploiting the results (16.15l) - (16.17l ) we arrive at the following formula for the 
derivative of the Neuman spectral zeta function 



(£")'(()) = (z")'(0) + Res£v(-i) 4 (In 2 - 1) f~ l (t)dt - 2 £ 



ln/(Q 
fit) 



dt 



I 2/+2 



+ ^4(0) + i^(0) [In f(a) + In /(ft)] - J ]r IfP^ (i + 1) #y, 2i+ i(a, ft) 

i=o y=o ^ v / 



£ 7, 2/+ 1 (a, ft) - J 7 ; 2i+ 1 (a, ft) 



-Res^(/ + i)| 2 ln2-g I A T 

IJ] 2i+l r s 

- YaYaY v ^ n (/) & 2 ' (a ' 6) + Res & (/) [^y+'-i^ 2i+1 (a ' ^ - 2 ' +1 (a ' f ' 

1=1 7=0 ^ > 



(6.18) 



where we have defined the functions 



Q jt ft) = Q u M, ft) + O] ^fi+Xa) + Q\ { (6)/^ +< (6) 



and 



Ty, ,-(a, ft) = Jy, ,-(a, 6) + <2y ffl)f 2i+i {a) In /(a) + 0+ t Q>)f 2i+i Q>) In /(ft) , 



(6.19) 



(6.20) 



with the assumption <27 +J .(a) - <3^ j .(ft) = 0. 

The anomalous scaling factor for Neuman boundary conditions can be obtained by referring, once again, 
to the relations (14. 14b . (14.151) and the explicit formula (16.141 ) for f (0). By noticing that Z^O) = 0, we obtain 

f"(0) - -Res&(-iJjT r 1 «df+i^(O)-|]Res^(0^#j li (fl,fe) . (6.21) 



VII. COEFFICIENTS OF THE HEAT KERNEL ASYMPTOTIC EXPANSION 

The analytic continuation of the spectral zeta function that we have found in section [TV] can be utilized 
in order to compute the coefficients of the heat kernel asymptotic expansion for the operator -A M on the 
manifold M. This follows from the fact that the spectral zeta function is related to the associated trace of 
the heat kernel through a Mellin transform. More specifically, one can prove that the following formulas 



hold (12,124 



31] 



^is-, = r(j)Resf(j), 



for s - n/2, (n - l)/2, • • • , 1/2 and s = -(21 + l)/2 for I e N , furthermore 

i-iy 



(7.1) 



(7.2) 
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for p e No- These relations are particularly useful because they produce any coefficient of the heat kernel 
asymptotic expansion in terms of either the residue or the value of the associated zeta function at a specific 
point. Since the manifold N is left unspecified the heat kernel coefficients on the warped product I Xf N 
will be expressed in terms of those on the manifold N. 

The explicit evaluation of the heat kernel coefficients associated with the operator -Am in arbitrary 
dimensions D is more conveniently performed by using the following formula 

(D-n\ (D-n\ 



which is valid for n < D. By keeping the dimension arbitrary, we are able to effectively compute all heat 
kernel coefficients which will then be given as a function of D. For Dirichlet boundary conditions, by 
choosing L - n — 1 in (14.21 ). Z{s) becomes an analytic function in the region %{s) > (D - n - l)/2 and 
therefore will not contribute to the value of the residue of £(s) at s = (D - n)/2. With this comment in mind, 
the only contributing factors are 

(D-n\ (D-n\ 1 (D - 1 - n\ ( D - 1 - n\ C b n „ , 

r (— ) ResA -' {—) - ut r Hr-r s( »Hr-)L f m • (1A) 

while from Aq(s) we have 

r(^)Re S A„(^) = -ir(^)Re S& (^)[/-» W + / D -»] , (7.5) 

and finally the contribution to the residue coming from A,(s) reads, when i = 2m + 1 with m e No, 

lD-n\ /D-n\ (D-n + 2m+l\ (D-n + 2m+Y 

r — — ResA 2m+1 — — = -T Res^V 



2 / \ 2 

:«/+: p/ Q-/i+2m+l l J\ 



V —, \ ; ' rf F j2m+l (t)f D - n ^ +2m+i (t)dt, (7.6) 

Z-i r(D-n+2m+l\ r ( m , ■ , 1\ I },-an+L\ jj w ' v ' 



and, when i = 2m with m e N + , 

(D-n\ (D-n\ (D-n \ (D-n \ 

r[—^—jResA 2m ^—^—] - -r\——+m)Res£ N \——+mj 

2m+l pj^n + j + m ) ( rh 



By recalling the relation (17.31 ) and by noticing that for the manifold /V this implies 

& N n=T — — Resfjv — — , (7.8) 
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(7.9) 



we obtain a general expression for the heat kernel coefficients on the manifold M when Dirichlet boundary 
conditions are imposed, namely 

# = f f D - n -\t)dt-\^_ l [f D ~"(a) + f D -'\b)] 

2 2 yjn 2 J a 4 — L J 

[t 2 ] 2i+2 Y (D-n+l + - + A fe 

6 2 1 ^r(^±i + /)r(/ + 7 - + i)X 

[t] 2/+1 p ( D^n + • + A , fe 



D-n+2 j+2i 



(a)} , 



where fflfL.n = for n = 0. The relation (I7.9I ) provides an efficient way of computing the coefficients of 
the heat kernel asymptotic expansion on the warped product I Xf N in terms of those of the manifold N. 
Furthermore it represents a very general result giving any heat kernel coefficient for an arbitrary dimension 
D. Obviously, more explicit results can be obtained once the manifold N and the warping function f(r) 
have been specified. It is instructive, at this point, to use ( T7.9I ) in order to provide the first few heat kernel 
coefficients, namely 



^ = f Am , 

2 J a 

[ f-\mt-\^[f{a)+f d ib)} , 

2 2 yn 2 J a 4 L J 



(7.10) 



(7.11) 



ff d -\t)d t -^ N l [f d ~\a) + f~\b)} + ^-^^ f f\i)f d -\t)dt , (7.12) 



and 



^n 3 



yn 2 J a 4 J 



2V^ 
1 

4 2 



f{b)f d - 2 {b)-f\a)f d -\a) 



•J a 



~ 16^ 



f"{b)f d -\b) + f'\a)f d ~\d) + ^-i (f'\b)f d -\b) + f'\a)f-\aj) 



(7.13) 



Higher order coefficients may be evaluated with the help of a simple computer program. We would like to 
point out that in order to obtain results in terms of geometric invariants of the manifold M in the expression 
for the heat kernel coefficients, a number of integration by parts in the variable t need to be performed. 
Also note that if N has a boundary, the above results clearly display corner contributions in the heat kernel 
coefficients. 
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The coefficients of the heat kernel asymptotic expansion for Neuman boundary conditions follow directly 
from those found in the Dirichlet case. In fact, by recalling (14- 14b . (14.151 ) and the formula (17.31 ) one obtains 



7 2 1 -T- L J 



2 

n-\ 



i y ( °- n+i + j) 

- Z hZ >w(o-U [QlMf D - n+2i+i M + ai i (b ) f D - n+2J+i (b)] , (7.i4) 

i=l j=0 1 [J+ 2) V {-^— ) 

which shows that the coefficients for the Neuman case can be expressed in terms of the coefficients associ- 
ated with the Dirichlet case. For completeness, the first four coefficients are provided below 



raNeu _ <2jDir 
^0 ~ ^0 ' 



j*f» = ^f + I< [f\a) +f d {b)} , 



(7.15) 



(7.16) 



and 



2 2 Z L J 2 Y^" 2 L J 



d-2 



f"{a)f d -\a) + f\b)f-\b) + — (/' 2 (a)/- 2 (a) + f\b)f d - 2 {b)) 



(7.17) 



VIII. CONCLUDING REMARKS 



In this work we have performed a detailed analysis of the analytic continuation of the spectral zeta func- 
tion associated with Laplace operators acting on scalar functions defined on the warped product manifold 
of the type I Xf N. We have exploited the explicit form of the analytic continuation in order to provide 
the zeta regularized functional determinant and the coefficients of the heat kernel asymptotic expansion as- 
sociated with the Laplacian on / Xy N. Since we have left the manifold ,/V unspecified the results for the 
functional determinant and the heat kernel coefficients have been given in terms of the spectral zeta function 
associated with the Laplace operator on the manifold 7Y. The technique for the analytic continuation that 
we have provided is based on the WKB asymptotic expansions of the eigenfunctions of the Laplacian on 
I Xf N. The novelty of this investigation resides in the generality of the technique presented which is valid 
for an arbitrary strictly positive warping function / and for any smooth compact cf -dimensional Riemannian 
manifold N. 

In this paper we have worked under the assumption that the warping function f{r) does not vanish at 
the endpoints of the interval /. It would be particularly interesting to generalize the analysis developed 
here in order to include the case of a vanishing warping function at either one or both the endpoints of 



26 



the interval I. We would like to point out that the results obtained in this work can be extended to cases 
in which I = [0, b] and f(r) ~ r as r — > 0. In this situation, the differential equation (12.81) has a regular 
singular point at r = and the WKB technique used to compute the uniform asymptotic expansion of the 
eigenfunctions can be directly exploited. The analysis proceeds in the same way as outlined in the previous 
sections by keeping the parameter a e (0, b). The analytic continuation of the spectral zeta function valid 
when I = [0, b] and f{r) ~ r as r -> is obtained by taking the limit a -> of (l4~2l and (l4~8l) - (14. 1 1 I) in the 
region of %(s) where the limit is well defined. Warped product manifolds for which f(r) ~ r as r — > are 
singular Riemannian manifolds which posses a conical singularity (like, for instance, the generalized cone 



and the spherical suspension tl4fl ). It would be very interesting to generalize this analysis to cases in which 
I = [0, b] and f(r) ~ r 6 when r — > with 6 + 1 as they would provide results for the spectral zeta function 
of Laplace operators on manifolds with singularities other than the conical one. 

Understanding the behavior of the Casimir energy for different geometric configurations has become 
a subject of major interest. It is well known that in the zeta function regularization scheme the Casimir 
energy is obtained by evaluating the zeta function associated with the problem under investigation at the 
point s = -1/2. Obviously this is done after a suitable analytic continuation of the spectral zeta function. 
The spectral zeta function that we have obtained in this paper could be utilized for the computation of 
the Casimir energy and force for a generalized piston configuration constructed from the warped product 
manifold / X/ N following the ideas developed for conical manifolds in IJj], 3]. In this case, however, a 
certain amount of numerical work would be necessary. The term Z(s) in (14.11) would contribute to the value 
of the spectral zeta function at the point s = -1/2 and it can only be evaluated numerically once a warping 
function has been specified. This types of investigations would help to shed some light on how the geometry 
influences the Casimir energy and force and it would be worthwhile to pursue research in that direction. 
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Appendix A: The functions 

In this appendix we list the functions F^jir) up to i = 3. By utilizing the recurrence relation (13.261) and 
the equality given in (13.271 ) we get 

d{d-2) f'\r) d(d-3)f(r)f"(r) d f\r) 
Fo,2ir) = — -y^ - y^- - -— , (Al) 
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{d 2 -2d- 12) /' 3 (r) (2rf + 9)f\r)f"(r) / (3) (r) 

Fi 2 (r) = 7 -. + — : — , (A2) 

%p(r) 8/ 4 (r) 8/ 3 (r) 

27/' 3 (r) 9f'(r)f"(r) 15 /' 3 (r) 

F 2, 2 W = -— y— + — — — , F 3j2 (r) = — , (A3) 

8/'(r) 8/ 6 (r) 8 / y (r) 

J(J 3 - Ad 2 - 20d + 48) y 4 (r) ^ 2 + 8 J - 24) / ,2 (r) f"(r) 
°' 3(r) " 128 / 4 (r) 32 T^W 

^-6)/" 2 (r) J(J-4)/'(r)/ (3) (r) f^(r) 

32 /2(r) 16 /2(r) 16 /(r) ' 1 



09^-38^-120) /' 4 (r) {\M 2 -54d- 144) f'\r)f"(r) 

p 1 = 

' 32 /6(r) 32 / 5 (r) 

3(rf + 3) f"\r) | (5rf + 12) f'(r)fW(r) 1 /< 4 >(r) 
16 / 4 (r) 16 /4(r) 16 /3(r) ' 



„ , . (25J 2 - 50d - 1014) f'\r) (25d + 366) f'\r)f"(r) 



64 f\r) 32 /7( r ) 

19 /" 2 (r) l f'{r)fV\r) 
32 /6(r) 8 f\r) ' 



(A6) 



663 /'V) 221 f\r)f'\r) 1105 /' 4 (r) 

^3^) = —^-— ^ , ^ 3 (r) = - — . (A7) 

Appendix B: The Functions £lk, i{a) 

By exploiting the cumulant expansion ( 13.311 ) we obtain 

M (a) = F M (a), * = 0,1,2 (Bl) 

flo,2(fl) = -^ F o, i(«) + ^o,3(«) , Oi,2(fl) - -Fo,i(a)Fi,i(a) + Fi f3 (a) , (B2) 

02.2(a) = j(a) - F , i(a)F 2 , ^a) + F 2>3 (a) , (B3) 

03.2(a) = -Fl, i(a)F 2 , i(a) + F 3>3 (a) , fi 4 , 2 (a) = -^F 2 j(a) + F 4>3 (a) , (B4) 

O 0>3 (a) - ^F 3 j(a) - F , i(a)F 0)3 (a) + F , 5 (a) , (B5) 

Oi, 3 (a) =F^ 1 (a)F 1)1 (a)-F 0i3 (a)Fi > i(a)-Fo,i(a)Fi i3 (a) + F li5 (a) , (B6) 
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"2,3(0) = ^o,i(a)^u(a)-^i,i(«)^i,3(a) + ^,i( a ) F 2,i(a) 



- F 0> 3 (a)F 2 , i (a) - F , i (a)F 2 , 3 (a) + F 2 , 5 (a) , 



(B7) 



^3,3(«) = -F\ A (a) + 2F QA (a)Fi A (a)F 2A (a) - Fi^(a)F 2A (a) 
- F h i(a)F 2>3 (a) - F 0> i(a)F 3>3 (a) + F 3>5 (a) , 



(B8) 



Q 4>3 (a) - Fj 1 (a)F 2 ,i(a) + F(j,i(a)F| 1 (a)-F2,i(a)F2,3(a) 



- F h i(a)F X3 (a) - F a i(a)F 4>3 (a) + F 4j5 (a) , 



(B9) 



Q 5;3 (a) = F u (a)F^ j(a) - F 2 , i(a)F 3>3 (a) - Fi, i(a)F 4>3 (a) + F 5j5 (a) , 



(BIO) 




F 2 j(a) - F 2> i(a)F 4> 3 (a) + F 6> 5 (a) . 
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